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Abstract 



We consider a model of Raman scattering for a two-dimensional S 



interaction. We observe a broadening of the line shape due to increased 



■ 1/2 Heisenberg Anti-Ferromagnet which includes a dynamical spin-phonon 

0^ 
oo 
O 

00 

Q ■ coupling with excited high-energy spin states. Our results are close to a model 

^ : 

0^ . of random, static exchange interactions, first introduced in this context by 



Haas et al. [J. Appl. Phys. 75, 6340, (1994)], which, when extended to large 



B 

. numbers of spins, explains experiments in the parent insulating compounds 

c : 

O ■ of high-Tc superconductors. 



I. INTRODUCTION 

Since the discovery of high-Tc superconductor Jil, the S = 1/2 quantum Heisenberg model 
has received considerable attention. This is largely due to well-accepted experimental ev- 
idence that suggests that these compounds can be described! by a two-dimensional (2D) 
doped Heisenberg model for spin S = 1/2. 

One of the experimental techniques used to study the excitations of these systems is Ra- 
man scattering. This technique was intensely used in the past for different antiferromagnets 
of spin S > 1. For such systems, the main features of the line shape were explained by 
Parkinson!, who used a spin-wave theory to account for the photon-two-magnon process 



1 



involved in the Raman scattering. Very good agreementcl was found for K2NiF4, which is 
well described by a spin S = 1 Heisenberg antiferro magnet. 

The results of Raman experiments in the parent insulating compounds of high-Tc su- 
perconductors, La2Cu04 and YBa2Cu306.2, show some qualitative differences with the line 
shape of previously studied antiferromagnetsi^. As in the case of K2NiF4, the line shape is 
centered at an energy corresponding to a spin-pair excitation. However, in contrast with the 
case of spin S = 1, the linewidth is very broad, and the spectrum has a very long tail that 
extends beyond the energy of four magnon excitations. Moreover, while the dominant con- 
tribution to scattering is in the so-called Big channel, there is also a significant contribution 
in the nominally forbidden Aig configuratiorfl. 

Since for lower spins the quantum fluctuations are larger, some theoretical analysis be- 
yond the mean-field spin-wave approximation has been attempted to explain the broad line 
shape. Numerical diagonalizations in a 4 x 4 lattice shows very little structure for the Big 
channel: essentially three peaks; a dominant two-magnon peak and two peaks identifiable as 
four-magnon excitationsi. This calculation gives vanishing line shape for the Aig channel. 
Although the structure of the line shape is clearly different from the one observed in the 
experiment, the first three moments of both lines are in good agreement. The first three mo- 
ments of the distribution have been calculated in good agreement to experiment by Singh et 
al. using cumulant expansions on a Heisenberg model with diagonal next-nearest-neighbor 
coupling^ Canali and Girvin0 used the Dyson-Maleev transformation taking into account 
processes of up to four magnons, and presented convincing evidence that the line shape 
cannot be explained by quantum fluctuations alone. Raman scattering Hamiltonians with 
a four-site exchange0 term have been proposed, these increase the relative weight of the 
four-magnon scattering, but does not explain the broadening of the peas. This is also con- 
sistent with the work of SugaiS, and Roger and DelrieuEi. Also, recent studies of spin-pair 
excitations in a spin 5* = 1 system, NiPSs, show a similar linewidth as those observed in the 
cuprate S0. 

From the above considerations, one concludes that it is necessary to invoke an additional 



process. It was emphasized by Nori et a/.El that the spin-phonon interaction can be respon- 
sible for the broad hnewidth, and for the finite cross-section in the otherwise absent Aig 
channel. They supported their arguments by computing the Raman cross section in finite 
lattices, using a Heisenberg model with random static exchange integrals, and averaging 
over configurations with equal weight. This calculation is in the spirit of the adiabatic or 
Born-Oppenheimer approximation^, that neglects the fluctuations of the phonon field. 

Adding phonons to the 2D S = 1/2 Heisenberg Hamiltonian (HAF) provides a 
mechanismlil for the otherwise forbidden Aig scattering, as well as allowing a coupling be- 
tween the ground and excited spin states. In a typical HAF, the phonon frequencyEl is about 
a third of the interaction constant J. The first excited state lies at 3J, so there is an order 
of magnitude energy difference between spin and phonon excitations. It was arguedEl that 
the separation of energy scales justifies the adiabatic phonon approximation. 

In the present work we consider the Raman scattering for a Heisenberg model with spin- 
phonon interaction, including the effects of the phonon djTiamics. We solve for the exact 
ground state of a system in which the vibrational degrees of freedom are allowed a finite 
number of excitations, an approximation valid for small but finite phonon frequencies. To 
the best of our knowledge, this is the first calculation of exact diagonalizations for a dynamic 
spin-phonon system. Some work including phonon dynamics beyond the adiabatic approxi- 
mation exists in the context of one-dimensional Peierls systems. Fradkin and Hirsh0 studied 
the electron-phonon interaction using quantum Monte Carlo simulations, and Proetto and 
Falicovli solved the case of two sites and one phonon. 

This paper is organized as follows. In Sec. II we present our Hamiltonian and scattering 
operator and discuss the rationale for our truncated-phonon model. In Sec. HI we present 
computational results and provide theoretical support for these results. We also present an 
alternative approach to including phonons and compare it to our model. Sec. IV is devoted 
to conclusions. 
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II. MODEL AND PROCEDURE 



We study a Heisenberg model with spin-phonon interaction described by the following 
Hamiltonian: 

H = {[j - o^ialj + (^ij) Si ■ Sj + ujoaljaij^ (1) 

<ij> 

where Si are spin 1/2 operators, creation (annihilation) operator for an Einstein 

phonon of frequency ujq, and the simbol "zj" refers to a link of a square lattice. In Jij = 
J — 6Jij we include the coupling between the phonons and the spin degrees of freedom 



through a "displacement" operator Xij = \J ^^^^'^Ij + parameter a in Eq.(|I|) is 

proportional to the spin-phonon coupling constant A = a ^ 2mujQ / h relating the change in 
the exchange integral 6Jij with the displacement: 6Jij = Xxij = a{alj + aij). 

Due to computational limitations, we need to restrict ourselves to a small number (6 and 
8) of spins. To have these as two dimensional as possible, we placed them in non-periodic 
ladder type structures as shown in Fig. 0. 

We use Einstein phonons to simplify the model, and we consider the highest-energy 
phonons as they are the closest in energy to the magnon excitations and will thus have the 
greatest coupling to the magnons. The occupation number of the phonon degrees of freedom 
at each link is in principle unrestricted; the number of phonons ranges from zero to infinity. 
This makes the problem intractable from the point of view of exact diagonalizations, since 
the resulting Hilbert space is infinite. In order to overcome this limitation, we restrict the 
possible number of phonons at each link by imposing the condition (a^)" = 0. Computational 
resources limit us to using n = 2,3. This approximation maps the phonon degree of freedom 
into a two- or three-state system at each link. The truncation of the phononic occupation 
states implies that the variations of Jij due to quantum fluctuations are bounded. For 
example, for n = 2, the maximum displacements for a given link is given by the "coherent" 
states \ip^) = 2~^/^{l±alj)\0) , and 6Jij = a{ip^\{alj + aij)\ilj^) = ±a. Quantum fluctuations 
themselves are certainly limited by our truncation, but since they are in general small for 
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small Uo, our approximation will account for the relevant dynamics in the regime SJij/J < 1 
and ojqI J < 1. 

We are interested in the Raman scattering intensity Ir{uj), given by: 

W = H mm\'Siuj-E, + E,), (2) 

V 

where \v) are the eigenstates and Ey the eigenvalues of H. We compute //j(c<j) by using the 
partial fraction expansion method of Gagliano and Balseiro0. The relevant operator R in 
Eq.(^ depends on the different configurations of the Raman experiment. In general!: 

-R = (-E'inc ■ Xij){Escat ' Xij)Si ■ Sj , (3) 

with Eiac and -Escat refer to the polarization of the incoming and scattered photon respec- 
tively, and Xij are unit vectors in the near-neighbor directions of the lattice. For square 
lattices, two common scattering configurations are the so-called Big and Aig symmetries. 
We let the lattice axis lie along the x and y directions, and we define x' and y' the directions 
along the lattice diagonals {x' = and y' = ^^), then Big corresponds to Escat \\ x' and 
Einc II y' 1 whereas Aig corresponds to both the incident and scattered photon polarized in 
the same direction: Einc \\ -Escat || a^'- 

For the case of the "pure" Heisenberg Hamiltonian (no phonons, or a = a;o = 0), the 
absorption corresponding to the Aig symmetry is zero at any non-zero frequency, since 
the Aig operator is proportional to the Hamiltonian. We show in the next section how 
the inclusion of the phonon dynamics gives rise to a finite absorption in this channel, and 
argue that a spin-phonon interaction accounts for the main features observed in the Raman 
experiments in the undoped copper oxides. Also, we compare the results obtained in the 
dynamical model with that of the Heisenberg model with random (static) exchange integrals, 
which was studied in this context by Nori et aH^ 
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III. RESULTS AND DISCUSSION 

In this section we present numerical results for systems of 6 and 8 spinsi§ corresponding 
to the geometry shown in Figure |l|. We first discuss the Aig symmetry. In this case, and 
in the absence of lattice distortions, the Raman operator commutes with the Heisenberg 
Hamiltonian. No line shape is observed in this case: Ir{uj) ~ 6(00). It was first pointed out 
by Nori et a/S that the presence of static disorder in the exchange integrals Jij changes the 
value of the commutator, and can give rise to a finite Aig signal quite similar to experiments. 
It was argued by Nori et a/.El that the spin-phonon coupling produces disorder in the values 
of Jij, in the limit where the vibrational motions are much slower than the spin degrees 
of freedom. Here, we are interested in how this limit is achieved in a system that — up to 
boundary effects — is translationally invariant, and includes the dynamics of the spin-phonon 
coupling. Consider our Hamiltonian Eq.([l|) and the operator R for the Aig symmetry. For 
c^o 7^ 0, [H, R] ^ and we expect a finite line shape. In Fig. ^ we show the I{uj) obtained 
for a = O.IJ and ujq = 0.05J. Some of the qualitative features of the experimental line 
shape are already present in these finite-size systems: there is a broad line shape of width 
~ 2J with a maximum at the two-magnon energy. Note that for our ladder geometries, 
since most sites have coordination number 3, the two-magnon energy in the Ising limit is 
located at 2 J, instead of the corresponding 3J of the square lattice. Due to the finite size of 
the lattice, the line shape consists of a series of peaks that are centered at the unperturbed 
energies, which are indicated by arrows in Fig. These unperturbed energies correspond 
to the manifold of two-and multi-magnon states. In turn, these "internal" peaks have a 
finite width that is due to the coupling with the phonons. The energy scale dominating this 
width is a. 

In order to test the validity of the Hilbert space truncation, we have increased the number 
of allowed excitations per phonon by one, with no qualitative changes in the line shape (see 
Fig. I) 

It is interesting to consider the following two limiting cases: (a) vanishing spin-phonon 
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coupling a ^ for finite plionon frequency uo, and (6) vanishingi phonon frequency ujq 
for finite a. In both limits the Aig line shape vanishes. In (a), the spins are uncoupled 
from the phonons, and the line shape vanishes because in this limit [H, R] = 0. Figure ^ 
illustrates how the line shape vanishes and shifts to lower energies as a — > 0. In limit (6), the 
phonons are "static", but remain coupled to the spins, because a ^ 0. Since the phonons 
have no dynamics, the eigenstates can be written in the form 

\^,{ujo = 0)) = \Mn})(^\{n}), (4) 

where |{^}) is a static configuration of displacements, in such a way that, for each link 
ij, one has (a|j + aij)\{n}) = ''7ij({'^})|{'^}), with f]ij{{n}) a c-number. The state \ipu{n}) 
involves spin degrees of freedom only, and is an eigenstate of the following Hamiltonian 

H{n} =T.[J- ( W)] ■ S, . (5) 

<ij> 

For ti^o = 0, the problem is that of an annealed configuration of displacements: for each 
configuration {n}, the couplings Jij are different, and one has to solve for the spin dynamics 
in the presence of this distribution of couplings. The complete spectrum is obtained by solv- 
ing Np decoupled eigenvalue problems, with Np being the phononic Hilbert space dimension. 
For our case of one excitation per link, Np = 2^""'^% with A^unks being the number of links 
of the lattice. What distinguishes this from a quenched disorder is the fact that the ground 
state of this problem corresponds to the lowest energy state of the spin system in an ordered 
background of couplings. Since this is also an eigenstate of R for the Aig symmetry, the line 
shape will be zero. In Fig. |^ we show how the line shape vanishes as uq decreases to zero. 
Note that from this analysis one concludes that, even in the case of cuq = 0, the line shape 
is non-zero at finite temperature. This is due to the presence of thermally excited states, 
which are eigenstates of a Hamiltonian that does not commute with R. The right column 
of Figure |^ shows the rescaled results, which illustrate the line shape approaching a limiting 
function of u whose overall amplitude vanishes as uq ^ 0. The main features of this limiting 
function should be accounted for by a model of static disordered couplings. In previous work 
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in the context of one-dimensional Peierls systemfll, it is argued that in the hmit of small 
phonon frequencies, the quantum lattice fluctuations can be modeled by a static, random 
Gaussian potential with zero mean. A model in which the exchange integrals are taken 
from a random configuration of couplings was also studied by Nori et a0. We can test this 
hypothesis in our dynamical model. Results of the comparison between the dynamical and a 
static disordered system are shown in Fig. Note that both the amplitude and the position 
of the peaks are in very good mutual agreement. The scattering peaks from the two models 
have their energies scaled with respect to each other. After they are rescaled, the location 
and magnitude of the peaks are very close for both models (see Fig. ^). The basis for this 
rescaling of the horizontal axis is the following. For our dynamic phonon model, the min- 
imum energy state corresponds to the lattice "maximally compressed": rjij = —{x)q. The 
tendency of the system to compress has the effect of renormalizing the interaction coefficient 
J, since the displacement Xij will tend to be a constant, and therefore J' = J + a{x)o- The 
average displacement will be greater as the possible number of excitations increases; this 
explains a shift in energy between 1 and 2 excitations as well. 

In the model with static disorder, the phonon energy goes to zero, but the fluctuations 
are still present. In the treatment of Nori et al^ the argument used is that the phonon 
energy is small compared to that of the spin excitations, and thus can be neglected. However, 
we have shown that for ujq strictly zero the line shape vanishes. At zero temperature, the 
disordered model should be compared with the dynamical model at infinitesimal Uq. We 
prove this statement by using perturbation theory in uq. Assume an ordering of the states 
for a = uq = 0, and let us label those states by the index z/. If a -C J and Uq = 0, each 
state \ipu) splits in a manifold I'ipuin}) of almost degenerate states of energies Eu{n}, that 
correspond to the "disordered" configurations. If we now turn on ujq, different states are 
going to mix, in such a way that the ground state can be written as 

iV'o) - iV'o Wo) + ^0 E c{M)\Mn}) ® \{n}) , (6) 

{n}^{n}o 

where {^^}o is the ordered configuration of couplings. A similar expansion can be used for 
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the excited states. If we use the fact that Aig R acts only on the spin degrees of freedom, 
and that different states \{n}) are orthogonal, to lowest order in uq, I{uj) will be given by 

JH ^ col E j:\iMn}\R\Mn})mco - E^n} + Eo{n}) . (7) 

{"} 

We have found computationally a falloff of I{uj) proportional to uJq. 

Now, if we keep the assumption of a <^ J, the states corresponding to each manifold 
are almost degenerate and we can approximate c{{n}) ~ (M — 1)^^^"^ ~ M^^^^, with M 
the number of configurations {n}. Therefore, in this limit I{uj) is given by an average over 
configurations with equal weight. 

We now present results for the Big configuration. In Fig. |^we compare the Big scattering 
from eight S = 1/2 Heisenberg spins for three different cases: bare HAF, with phonon 
interactions, and with static disorder. 

We have argued that peaks at higher excitations will eventually merge into one broad 
peak as the lattice size increases. To show that this is the case, we have calculated the 
Aig and Big peaks for a 3 x 4 non-periodic lattice with static disorder. The results are 
shown in Fig. |^. These are in good qualitative agreement with experiments, except that 
the experimental scattering is centered at a; ~ 3 J instead of 2J in our case. This is due to 
finite size effects and to the lattice being non-periodic. 

The model with static disorder, when extended to larger numbers of spins than can be 
obtained with the dynamic phonon model, begins to duplicate experiments. The agreement 
between the static and dynamic phonon models suggests that the dynamic phonon model 
will also agree with experiments. 

IV. CONCLUSIONS 

In this paper we have presented a model of phonon-magnon interaction with a truncated 
phonon space. We discussed how this model can explain the otherwise forbidden Aig scat- 
tering and showed that it does in fact give rise to Aig scattering. We discussed how coupling 
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to spin excitations leads to a broadening of the line shape. We have proven that phonons can 
be modeled by static disorder and compared the results obtained by using quenched phonons 
to our dynamical model. We showed that the amplitude of the Aig scattering should fall off 
as (Vq. We showed that both the Aig and Big scattering is broadened. In conclusion, the 
broadening is due to coupling to spin excitation states through the phonon interaction. 
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FIGURES 

FIG. 1. Ladder geometries used in our numerical computations. Each dot represents 
a spin 1/2. Each sohd hne connects nearest-neighboring spins, and represents a Heisenberg 
Anti-Ferromagnetic couphng modulated by an Einstein phonon. There is one phonon per cou- 
pling. 

FIG. 2. Raman cross section I{uj) versus u/J for Aig scattering. The parameters of the 
dynamical Hamiltonian are a = O.IJ, loq = 0.05 J, and the number of phonon excitations per link, 
n, is 2. The vertical scale is arbitrary. The arrows are located at the eigenenergies of the model 
without phonons corresponding to eigenstates of total spin S = 0. 

FIG. 3. Comparison between one-and two-phonon (lighter curve) excitations per link, for the 
Raman line shape in the Aig symmetry. The parameters are a = O.IJ, and ujq = 0.05J, for 
a ladder of 6 spins. 

FIG. 4. Raman line shape for the Aig symmetry. Here we consider 6 spins (including 1 
excitation per link), and ujq = 0.05 J. The strength of the spin-phonon coupling is varied: (a) 
a = 0.2J, (b) a = O.IJ (c) a = 0.05J, (d) a = 0.025J. The vertical scale is arbitrary but is the 
same for the different values of a. 

FIG. 5. Raman spectrum versus uj/J for the Aig symmetry. These results correspond to 6 
spins coupled to one phonon excitation and a = O.IJ. (a)-(c) are all displayed on the same scale, 
(d)-(f) have been rescaled. ujq = 0.3 for (a) and (d), = 0.1 for (b) and (e), and ujo = 0.025 for 
(c) and (f). 

FIG. 6. Raman spectra for the Aig symmetry for 8 spins coupled to one phonon excitation 
and also scattering from an average over configurations of disorder, (a) a = O.IJ, loq = 0.05J, and 
disorder chosen from a Gaussian with variance a = 0.2J. (b) a = 0.3J, uq = 0.3J, and disorder 
chosen from a Gaussian with variance a = 0.5J. 
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FIG. 7. Raman scattering for the Big symmetry using (a) the bare Heisenberg Hamiltonian 
(a = 0), (b) with dynamical phonons a = O.IJ, = 0.05J, (c) with dynamical phonons a = 0.3J, 
loq = 0.3J. The two bottom curves correspond to static Gaussian disorder of variance a = 0.2J 
(d), and a = 0.5J (e). 

FIG. 8. Raman scattering from a 3 x 4 non-periodic lattice in the Aig and Big geometries. 
The spin interactions have static Gaussian disorder of variance a = 0.5 
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